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In this paper we extend the color dipole formalism for the study of leading neutron production 
ine+p—s-e + n + X collisions at high energies and estimate the related observables which were 
measured at HERA and could be analysed in future electron-proton (ep) colliders. In particular, we 
calculate the Feynman xf distribution of leading neutrons, which is expressed in terms of the pion 
flux and the photon-pion total cross section. In the color dipole formalism, the photon-pion cross 
section is described in terms of the dipole-pion scattering amplitude, which contains information 
about the QCD dynamics at high energies and gluon saturation effects. We consider different 
models for the scattering amplitude, which have been used to describe the inclusive and diffractive 
ep HERA data. Moreover, the model dependence of our predictions with the description of the pion 
flux is analysed in detail. We demonstrate the recently released HI leading neutron spectra can be 
described using the color dipole formalism and that these spectra could help us to observe more 
clearly gluon saturation effects in future ep colliders. 

PACS numbers: 12.38.-t, 24.85.-|-p, 25.30.-c 
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I. INTRODUCTION 

In high energy collisions the outgoing baryons which have large fractional longitudinal momentum (xp > 0.3) and 
the same valence quarks (or at least one of them) as the incoming particles are called leading particles (LP). The 
momentum spectra of leading particles have been measured already some time ago [l[ . Very recently, high precision 
data on leading neutrons produced in electron-proton reactions at HERA at high energies became available 0 . 

Leading neutron production is a very interesting process. In spite of more than ten years of intense experimental 
and theoretical efforts, the xp (Feynman momentum) distribution of leading neutrons remains without a satisfactory 
theoretical description [3l-[l^. Monte Carlo studies, using standard deep inelastic scattering (DIS) generators show 
SEl that these processes have a rate of neutron production a factor of three lower than the data and produce a 
neutron energy spectrum with the wrong shape, peaking at values of xp below 0.3. In order to fit the data the existing 
models need to combine different ingredients including pion exchange, reggeon exchange, baryon resonance excitation 
and decay and independent fragmentation [Ml. Moreover, the incoming photon interacts with the pion emitted 
by proton and then rescatters, interacting also with the emerging neutron and giving origin to significant absorptive 
corrections, which are difficult to calculate 0[ll. As it can be seen in Fig. [U this process is essentially composed 
by a soft pion (or reggeon) emission and by the subsequent photon-pion interaction at high energies. Pion emission 
has been studied for a long time and according to the traditional wisdom it can be described by a simple interaction 
Lagrangian of the form gijj'y^'K'ip, where if and tt are the nucleon and pion field respectively. The corresponding pion- 
nucleon amplitude must be supplemented with a form factor, which represents the extended nature of hadrons and 
at the same time regularizes divergent integrals. The precise functional form of the form factor is ver y h ard (if not 
impossible) to extract from first principle calculations. We must then resort to models. Very recently it has 

been argued that the Lagrangian which is more consistent with the chiral symmetry requirements is the one with a 
pseudovector couplings. One of the goals of the present study is to test in a new context the better founded splitting 
function, f{y), derived in M- Additionally, forward hadron production is very important also for cosmic ray physics. 
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FIG. 1: Leading neutron n production in ep enX interactions at high energies, xl is momentum fraction of the proton 
carried by the neutron n. 


where highly energetic protons reach the top of the atmosphere and undergo successive high energy scatterings off light 
nuclei in the air. In each of these collisions, a projectile proton (the leading baryon) looses energy, creating showers 
of particles, and goes to the next scattering. The interpretation of cosmic data depends on the accurate knowledge 
of the leading baryon momentum spectrum and its energy dependence. A crucial question of practical importance 
is the existence or non-existence of the Feynman scaling, which says that the xp spectra of secondaries are energy 
independent. In cosmic ray applications we are sensitive essentially to the large xp region (the fragmentation region), 
which probes the low Bjorken-cc component of the target wave function. In this kinematical range nonlinear effects 
are expected to be present in the description of the QCD dynamics (for recent reviews see [l^), associated to the 
high parton density. The state-of-art framework to treat QCD at high energies is the Color Glass Condensate (CGC) 
formalism [13, which predicts gluon saturation at small-x, with the evolution with the energy being described by an 
infinite hierarchy of coupled equations for the correlators of Wilson lines - the Balitsky-JIMWLK hierarchy [T^. In 
the mean field approximation, this set of equations can be approximated by the Balitsky-Kovchegov (BK) equation 
Cliil, which describes the evolution of the dipole-target scattering amplitude with the rapidity Y = ln(I/a;). 

In this paper we propose to treat the leading neutron production at HERA using the color dipole formalism, which 
is able to describe the inclusive and diffractive HERA ep data. Our goal is to extend this successful formalism, which 
have its main parameters well determined, for leading neutron production. As a consequence, our predictions are 
free parameter, depending only from the choices for the models of the pion flux and dipole scattering amplitude. 
Moreover, the use of the color dipole formalism allows to estimate the contribution of the gluon saturation effects for 
the leading neutron production in the kinematical range which was probed by HERA and which will probed in future 
electron-proton colliders. Einally, it allows to investigate the relation between the Eeynman scaling (or its violation) 
and the description of the QGD dynamics at high energies. It is important to emphasize that, in the near future, 
Eeynman scaling will be investigated experimentally at the LHG by the LHGf Collaboration [20l4^ . 

This paper is organized as follows. In the next Section we present a brief review of the leading neutron production 
in ep collisions as well as the different models for the pion flux are discussed. Moreover, the treatment of the process 
using the color dipole formalism is presented and the main assumptions are analysed. In Section IlIIl we analyse the 
dependence of our predictions in the pion flux and in the scattering amplitude. A comparison with the recent HI 
data is presented and the Feynman scaling is analysed. Finally, in Section llVI we summarize our main conclusions. 


II. LEADING NEUTRON PRODUCTION IN THE COLOR DIPOLE FORMALISM 

A. The cross section 


Let us review the main formulas of leading neutron production. At high energies, this scattering can be seen as 
a set of three factorizable subprocesses (See Fig. [IJ: i) the photon fluctuates into a quark-antiquark pair (the color 
dipole), ii) the color dipole interacts with the pion, present in the wave function of the incident proton, and iii) the 
leading neutron is formed. The differential cross section reads: 


d'^a{W,Q'^,XL,t) 




dxpdt 


( 1 ) 
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where is the virtuality of the exchanged photon, W is the center-of-mass energy of the virtual photon-pion system. 
It can be written as = {1 — Xl) W'^, where W is the center-of-mass energy of the virtual photon-proton system. 
As it can be seen in Fig. [U is the proton momentum fraction carried by the neutron and t is the square of the 
four-momentum of the exchanged pion. In terms of the measured quantities xl and transverse momentum pT, the 
pion virtuality is: 

^_ Pt (1 ~ ^L){ml - mlxL) ^ 2 ) 

XL XL 

The flux of virtual pions emitted by the proton is represented by and Q"^) is the cross section for the 

interaction between the virtual-photon and the virtual-pion at center-of-mass energy W. 


B. The pion flux 


The pion flux fTr/p{xL,t) (also called sometimes pion splitting function) is the virtual pion momentum distribution 
in a physical nucleon (the bare nucleon plus the “pion cloud”). It was first calculated in the early studies [2^ of deep 
inelastic scattering (DIS), where a pseudoscalar nucleoii-pion-nucleon vertex was added to the standard DIS diagram. 
These early calculations were further refined in Refs. and even extended to the strange and charm sector [^ . 
Since pion emission is a nonperturbative process, once we depart from the single nucleon state, we should consider 
a whole tower of meson-baryon states, having to deal with a series for which there is no rigorous truncation scheme. 
This led some authors to use light cone models [2^ for the pion momentum distribution, where the dynamical origin 
of the pion is not mentioned and one tries to determine phenomenologically the relative weight of the higher Fock 
states. 

In all the calculations of the pion flux a form factor was introduced to represent the non pointlike nature of hadrons 
and hadronic vertices, which contain a cut-off parameter determined by fitting data. The most frequently used 
parametrizations of the pion flux (Mi have the following general form: 

U/pixL,t) = - XL)^-^‘^^^^[F{xL,t)]'^ (3) 


where gp^p/{‘iTr) = 14.4 is the ir^pp coupling constant, TO,r is the pion mass and a{t) will be defined below. The form 
factor F{xL,t) accounts for the finite size of the nucleon and pion. We will consider the following parametrizations 
of the form factor: 


Fi{xL,t) = exp 


2 {t-mlY 
. a-XL). 


a{t) = 0 


( 4 ) 


from Ref. where i? = 0.6 GeV ^. 


F2{xL,t) = 1 , a{t) = a(t)^ (5) 

from Ref. Q, where a,r(t) — t (with t in GeV^) is the Regge trajectory of the pion. 

F3{xL,t) = eTq)[b{t - ml)] , a{t) = a{t)^ ( 6 ) 

from Ref. Q, where b = 0.3 GeV“^. 

F,{xL,t) = , a{t )=0 (7) 

from Ref. Q, where = 0.74 GeV. 


F5{xL,t) 




a{t) = 0 


( 8 ) 


also from Ref. Q, where A^^ = 1.2 GeV. In the case of the more familiar exponential monopole © and dipole 
(|8|) forms factors, the cut-off parameters have been determined by fitting low energy data on nucleon and nuclear 
reactions and also data on deep inelastic scattering and structure functions. 
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In all analyses of pion cloud effects in DIS the authors have used pseudoscalar (PS) coupling, which is, by itself, 
inconsistent with chiral symmetry. More recently, in the authors used an effective chiral Lagrangian for the 
interaction of pions and nucleons consistent with chiral symmetry. Unlike the earlier chiral effective theory calculations 
which only computed the light-cone distributions of pions or considered the non analytic behaviour of their moments 
to lowest order in the pion mass, in the authors have computed the complete set of diagrams relevant for DIS from 
nucleons dressed by pions resulting from their Lagrangians, without taking the heavy baryon limit. In particular, 
they have demonstrated explicitly the consistency of the computed distribution functions with electromagnetic gauge 
invariance. In [Tsil they applied the previously computed pion distribution to the study of the d — u asymmetry in 
the nucleon. For our purposes, the relevant contribution of the pion momentum distribution reads M- 


f-ir/piy) 


9Arnl 

(47r/^)2 



, 2 y{PT + 


(9) 


where y = 1 — xl, Pa = 1.267 is the nucleon axial charge, /tt = 93 MeV and Ac = 0.2 GeV. In order to obtain the 
final cross section we multiply (O by 2 (see 0[ll) and insert it into © after integrating the latter over t (or over 
Pt)- 

In what follows, we shall use the pion fluxes listed above, in equations ©-([9]), denoting them by /i, / 2 , ... /e, 
respectively. The main purpose of our calculation will be to show that the dipole approach gives a good description of 
data and the pion flux is just one element of the calculation. Therefore we will make no effort to choose one particular 
form. Nevertheless, we note that one important phenomenological constraint that these pion fluxes must satisfy is to 
reproduce the d — u asymmetry in the proton sea measured by the E866 Collaboration [T7I | . Among the pion fluxes 
mentioned above only ([5]) (see [^) and © have been confronted with these data. As it will be seen later, the results 
discussed here are very sensitive to the choice of the pion flux and hence leading neutron spectra may be used to 
constrain its shape. 

Before closing this subsection we would like to mention that the diagram shown in Fig|T]represents only the dominant 
contribution to leading neutron production. Other isovector meson exchanges, such as p or 02 , can also contribute to 
the leading neutron spectrum. Moreover the p —>■ A transition can also contribute to neutron production through the 
subsequent decay A —>• mr. Theoretical studies show that processes other than direct pion exchange are expected to 
give a contribution of less than 25 % of the cross section (3-0, HI] . The effect of these backgrounds to the one pion 
exchange is to increase the rate of neutron production. However this effect is partially compensated by the absorptive 
rescattering of the neutron, which decreases the neutron rate by approximately the same amount [IlL [l2| . 


C. The photon-pion cross section 

In order to obtain the photon-pion cross section we will use the color dipole formalism, as usually done in high 
energy deep inelastic scattering off a nucleon target. In this formalism, the cross section is factorized in terms of the 
photon wave functions 4', which describes the photon splitting in a qq pair, and the dipole-pion cross section adi^- We 
have 

[ dz f d'^r'^\'i>T,L{z,r,Q'^)\^ adT,{x,r) ( 10 ) 

«/ 0 ^ T T' 


where 

-I- TOj -I- 

^ ~ -t Q2 “ (I - xl)W'^ + Q2 


( 11 ) 


is the scaled Bjorken variable and the variable r defines the relative transverse separation of the pair (dipole). 
Moreover, we have that the photon wave functions are given by 


/ 


( 12 ) 

(13) 


for a longitudinally (L) and transversely (T) polarized photon, respectively. In the above expressions = z{\ — 
z)Q'^ + TOj , Kq and Ki are modified Bessel functions and the sum is over quarks of flavour / with a corresponding 
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quark mass m/. As usual z stands for the longitudinal photon momentum fraction carried by the quark and 1 — z is 
the longitudinal photon momentum fraction of the antiquark. 

The main input in the calculations of is the dipole-pion cross section. In what follows, for simplicity, we will 
assume the validity of the additive quark model, which allows us to relate Od-K with the dipole-proton cross section, 
usually probed in the typical inclusive and exclusive processes at HERA. Basically, we will assume that 


OdT,{x,r) = - • cjdp{x,r) 


(14) 


This is assumption is supported by the study of the pion structure function in the low x re gime presented in . It 
also gives a good description of the previous ZEUS leading neutron spectra, as shown in [a, [l3- On the other hand, 
the direct application of o to HERA photoproduction data [3l| leads to the result /o^ = 0.32 ± 0.03, which 
is factor 2 lower than the ratio given above. This subject certainly deserves more detailed studies in the future. For 
our purposes, the use of relation (I14p allows us to estimate Ud-n without additional free parameters. In the eikonal 
approximation the dipole-proton cross section adp is given by: 


adp{x,r) = 2 


cPbM^{x, r, b) , 


(15) 


where r, b) is the imaginary part of the forward amplitude for the scattering between a small dipole (a colorless 

quark-antiquark pair) and a dense hadron target, at a given rapidity interval Y = ln(l/a;). The dipole has transverse 
size given by the vector r = x — y, where x and y are the transverse vectors for the quark and antiquark, respectively, 
and impact parameter b = {x + y)j2. As mentioned in the introduction, at very high energies (and very low x) the 
evolution with the rapidity Y of ^/^{x, r, b) is given by the Balitsky-Kovchegov (BK) equation assuming the 

translational invariance approximation, which implies J\fP{x, r, b) = J\fP{x, r)S{b) and adp{x, r) = ao ■ M^{x, r), with 
the normalization of the dipole cross section (ctq) being fitted to data. Alternatively, we can describe the scattering 
amplitude J\fP{x,r) using phenomenological models based on saturation physics constructed taking into account the 
analytical solutions of the BK equation which are known in the low and high density regimes. The main advantage 
in the use of phenomenological models is that we can easily to compare the linear and nonlinear predictions, which is 
useful to determine the contribution of the saturation effects for the process under analysis. In what follows we will 
consider as input the phenomenological models proposed in Refs. 133 - 1^ . In particular, the BMS model, proposed 
in Ref. [s^ and updated in [s^, was constructed so as to reproduce two limits of the LO BK equation analytically 
under control: the solution of the BFKL equation for small dipole sizes, r <C 1/Qs{x)^ and the Levin-Tuchin law for 
larger ones, r ^ \/Qs{x). In the updated version of this parametrization [s^, the free parameters were obtained by 
fitting the new HI and ZEUS data. In this parametrization the forward dipole-proton scattering amplitude is given 
by 


N^{x, r) 



for rQs {x) <2 , 
for rQs{x) > 2 , 


(16) 


where a and b are determined by continuity conditions at rQs(x) = 2, 7 s = 0.7376, k = 9.9, Afo = 0.7 and Qs is the 
saturation scale given by: 


Qi(^)=Qi 


( 17 ) 


with xq = 1.632 X 10 A = 0.2197, = 1.0 GeV^. The first line of Eq. (ITHll describes the linear regime whereas 

the second one includes saturation effects. On the other hand, in the GBW model [Hj, the dipole-proton scattering 
amplitude is given by 


Af^{x, r) = 1 — exp 



(18) 


with Qq = 1.0 GeV^, Xq = 3 x 10 ^ and A = 0.288. Finally, we will also use the dipole-proton cross section estimated 
with the help of the DGLAP analysis of the gluon distribution, which is given by [s^: 

TT^ 

cTdipix, r) = —r'^asxg{x, lO/r^) (19) 

where xg{x,Q^) is the target gluon distribution, for which we use the GTEQ6 parametrization [s^. The above 
expression represents the linear regime of QCD and is a baseline for comparison with the nonlinear predictions. 
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FIG. 2: Scaled Bjorken variable x, (Eq. Hill , as a function oi xl for two different virtualities: (a) — 10 GeV^ and (b) 

= 53 GeV^ 


III. RESULTS 


As it was seen in the previous sections, the leading neutron spectrum depends essentially on two main ingredients: 
the pion flux and the dipole-pion cross section. The latter is very sensitive to the value of the involved Bjorken x, or, 
in our case, x. In Fig. [5] we show a; as a function of Xl- From the figure we can conclude that at the highest values 
of the present photon-hadron energies {W) and at the lowest values of we enter deeply in the low x domain. This 
will be even more so if measurements can be carried out at higher energies, but already at the available energies we 
can see that the leading neutron spectrum receives a contribution from the kinematical range where the nonlinear 
effects are expected to be present and hence it is interesting to investigate their influence. 

We address now the dependence of our predictions for the leading neutron spectra on the models of the pion flux 
and of the dipole-pion scattering amplitude. It is important to emphasize that the spectrum is proportional to the 
product of these two quantities. Consequently, in what follows we will initially assume a given model for one of these 
quantities and analyse the dependence on the other quantity. Having done that, we will choose one combination, 
calculate the normalized cross section (l/aDis)dcr/dxL and compare with the experimental data. The normalization 
was taken from Ref. [s^ and is given by: 


ctdis = 


An^aem c 
xd 


( 20 ) 


where c = 0.18, /3 = d ln((3^/Ag) with d = 0.0481 and Aq = 0.292 GeV. In Fig. [3] (a) we consider the IIMS 
model for the scattering amplitude and estimate the leading neutron spectrum using the different models of the pion 
flux discussed in the previous section. We observe that the behavior at medium and small values of xl is strongly 
dependent on the choice of the model. Using the cut-off Ac in the range defined in Ref. [l^, our results suggest 
that model fe is disfavored. In Fig. [3] (b) we consider the /2 model for the pion flux and estimate the spectra using 
different models of the scattering amplitude. We can see that the predictions have similar behavior in the transition 
between the small and large-a;/, regimes, but with the magnitude being dependent on the description of the QCD 
dynamics. In the previous analysis we have considered = 53 GeV^. However, the experimental data were obtained 
for a wide range of virtualities, which implies that there is some arbitrariness present in the choice of the value of 
used in our calculations. In order to estimate the uncertainty present in this choice, in Fig. [3] (c) we consider the 
/2 and IIMS models and calculate the spectrum for different values of the photon virtuality. We observe that our 
predictions are more compatible with the experimental data if larger values of are assumed. 

So far we have assumed that the photon hits only the pion, in a type of “impulse aproximation”. However it has 
been shown in [l^ that very often the photon hits also the neutron, specially in the low domain, where the photon 
has low resolving power. In these cases the extra interactions generate the so called absor ptiv e corrections, which 
can be estimated with models. In the case of leading neutron production, it has been show n 11211 that the corrections 
are not very large and affect almost uniformly all the xl spectrum. Using the results of |12j |. we shall here simply 
multiply our spectra by the absorptive correction factor, Kabs =0.7. 

In Fig. |4]we give our description of data, which is quite reasonable. All the parameters contained in the dipole 
cross section have been already fixed by the analysis of other DIS data from HERA. The pion flux has been fixed 
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(c) 


FIG. 3: Leading neutron spectra compared with experimental data a.t W = 100 GeV considering (a) the IIMS scattering 
amplitude and different choices of the pion flux = 53 GeV^); (b) the /2 pion flux and different choices for the scattering 
amplitude (Q^ = 53 GeV^); and (c) the /2 pion flux and the IIMS scattering amplitude for different values of photon virtuality 
































FIG. 4: Leading neutron xl spectrum for three photon-proton energies compared with recent HI data. 


by fitting data of low energy hadronic collisions. As metioned above, it has been argued that the factor | predicted 
by the additive quark model to relate the photon-pion and photon-proton cross sections could be smaller, closer to 
i. Using this different factor, we would then underpredict the data in the low xl region. In principle, this would 
be acceptable since in this region other processes are expected to play a role. As already emphasized, this subject 
deserves more detailed studies. 

Finally, we analyse the Feynman scaling in the leading neutron spectra and the contribution of nonlinear effects for 
this process. As mentioned before this process is surprisingly sensitive to low x physics. In particular, the large Xl 
{xl > 0.3) we observe the transition from the large to small x domain. It is therefore interesting to check whether 
gluon saturation effects are already playing a significant role. As it can be seen from Eqs. H]) and (fTUl) . all the energy 
dependence is contained in the dipole forward scattering amplitude M. Therefore, we can estimate the contribution 
of the nonlinear effects comparing the results obtained with saturation model, e.g. the IIMS model given in (II6L 
with those obtained with the linear model, Eq. (HU). Moreover, the theoretical expectation can be obtained using 
the GBW model for the scattering amplitude, Eq. m- In the linear limit, when the dipole radius is very small (or 
equivalently is very large) or the saturation scale is very small (and hence the energy is not very high), we can 
expand the exponent and obtain 


crd7v{r, x) oc (To A/'(r, x) ~ ag 


Q^jxy 


~ croQ^o^ 


(I-xl)IU2+Q2 


Q2 _i_ 


f 


( 21 ) 


Consequently, in this regime we see that the leading neutron xl spectrum will depend on W. In a complementary 
way, in the nonlinear limit, when the dipole radius is very large (or equivalently is very small) or the saturation 
scale is very large (and hence the energy is very high), we obtain 

adTv (r, x) oc (To A/'(r, x) ^ ag (22) 

which is energy independent. One could argue that the mere inspection of Eq. (HD would suggest that at some asymp¬ 
totically high energy the photon-pion cross section would reach some “black disk” limit and the energy dependence 
would disappear. We would like to emphasize that the information contained in (EH) and (1221) is much richer and 
indicates the route through which the asymptotic limit is reached and the role played by nonlinear effects. These 
expectations can be compared with those obtained using the IIMS and DGLAP models for the dipole-pion cross 
section. In Eig. [5] (a) we show the spectra obtained in a purely linear approach. As expected we see a noticeable 
energy dependence. In contrast, the nonlinear predictions presented in Eig. [5] (b) show a remarkable suppression of 
the energy dependence at low values of consistent with the expectations. These results indicate that the Feynman 
scaling (and how it is violated) can be directly related to the QCD dynamics at small-a:. 
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(a) (b) 


FIG. 5: (a) Leading neutron spectra for different energies considering (a) the linear DGLAP model and (b) the nonlinear IIMS 
model for the dipole-pion cross section. 


IV. SUMMARY 

In this work we have studied leading neutron production in e + p —> e + n X collisions at high energies and 
have calculated the Feynman xl distribution of these neutrons. The differential cross section was written in terms 
of the pion flux and the photon-pion total cross section. We have proposed to describe this process using the color 
dipole formalism and, assuming the validity of the additive quark model, we have related the dipole-pion with the 
well determined dipole-proton cross section. In this formalism we have been able to estimate the dependence of the 
predictions on the description of the QCD dynamics at high energies as well as the contribution of the gluon saturation 
effects for the leading neutron production. With the parameters strongly constrained by other phenomenological 
information, we were able to reproduce the recently released HI leading neutron spectra. One of our most interesting 
conclusions is that leading neutron spectra can be used to probe the low x content of the pion target and hence it 
is a new observable where we can look for gluon saturation effects. At higher energies this statement will become 
more valid. Moreover saturation physics provides a precise route to Feynman scaling, which will eventually happen 
at higher energies. These conclusions motivate more detailed studies, in particular in the subjects which are the main 
source of theoretical uncertainties in our calculations: i) the validity of the additive quark model; ii) the factor Kabs 
accounting for the absorptive corrections, which are model dependent; iii) the precise form of the pion flux; iv) the 
precise form of the dipole cross section. It is important to emphasize that none of these numbers or functions is 
free. On the contrary, they are subject to severe constraints from other experimental information and the freedom to 
choose them will be further reduced in the future. 
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